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ABSTRACT. One might at first tend to think that the growth of plants and animals, because 
of their elaborate forms, are ruled by highly complex laws. However, this is surprisingly 
not always true: many aspects of the growth of plants and animals may be described by 
remarkably simple mathematical laws. An obvious example of this are the seashells and 
snails, as we show here: with a very simple model it is possible to describe and generate 
any of the many types of seashells of the classes of Gastropods, Bivalves, Cephalopods and 
Scaphopods that one may find classified in a good seashell bookguide. 


«Beauty of style and harmony and grace and good rhythm 
depends on simplicity.> 


— PLATO 


<There is much beauty in nature’s clues, and we can all rec- 
ognize it without any mathematical training. There is beauty, 
too, in the mathematical stories that start from the clues and 
deduce the underlying rules and regularities, but it is a dif- 
ferent kind of beauty, applying to ideas rather than things. 
Mathematics is to nature as Sherlock Holmes is to evidence.> 

— I. STEWART 


1. How SEASHELLS GROW 


The idea that mathematics is deeply implied in the natural forms goes back to the Ancient 
Greeks. Many aspects of the growth of plants and animals may be described by remarkably 
simple mathematical laws, in spite of their elaborate forms (cf., for instance, the classical 
book of D’Arcy Thompson [9] and the recent book of Stephen Wolfram [10]). 

An obvious example of this are the seashells and snails [6]. Why do so many shells form 
spirals? As far as the animal that lives in a shell grows he needs the shell to grow with him 
in order to keep accommodating him. The fact that the animal which lives at the open edge 
of the shell places new shell material always in that edge, and faster on one side than the 
other, makes the shell to grow in a spiral. The rates at which shell material is secreted at 
different points of the open edge are presumably determined by the anatomy of the animal. 
And, surprisingly, even fairly small changes in such rates can have quite tremendous effects 
on the overall shape of the shell, which is in the origin of the existence of a great diversity 
of shells. 
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A two-dimensional version of this fact may be observed in the growth of horns. Like nails 
and hair, a horn grows by the addition of material as its base. In order to get a perfectly 
straight structure, the rate at which material is added must be exactly the same on each 
side of the base: 


On the other hand, if there is some difference (indicated in percentage, in the figures 
below), one of the edges of the horn will be longer than the other and, inevitably, the horn 
will coil to the side where less material is added, following a spiral curve: 


wow £ 


10% 100% 


It is essentially a three-dimensional version of this phenomena that yields the spiral 
structures of the shells of mollusks. Besides that, the mollusk does not enlarge its shell in 
a uniform way: it only adds material in one of the edges of the shell (the open or “growth” 
ending) and makes it in such a way that the new shell is always an exact model, to scale, 
of the smaller shell. This growth process yields an elegant spiral structure (very visible 
when the shell is sliced). The widths of the straight lines that link the shell center (the 
spiral origin) to the points of the shell increase, but the amplitudes of the angles defined by 
those lines and the corresponding tangents to the shell are constant, that is, shells follow 


an equiangular spiral: 


This fact was identified in the 17th century by Christopher Wren. A clear mathematical 
model of shell growth modes based on equiangular spirals was given by Henry Moseley 
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in 1838, and the model used here is a direct extension of his (M. B. Cortie [I]). Careful 
studies from the mid-1800s to mid-1900s validated Moseley’s basic model for a wide variety 
of shells. 

As indicated in the figure above, given a point O (the origin or the pole of the spiral), 
an equiangular spiral is a curve such that the angle a (the equiangular angle) between the 
tangent in each point P of the curve and the radial line OP is constant. Jacob Bernoulli 
(1654-1705) called it the Spira mirabilis (the marvellous spiral). It was first described 
mathematically in 1638 by René Descartes (who believed that “only mathematics is certain, 
so everything must be based on mathematics”). Its parametric equation] is given, in polar 
coordinates r and 0, by 

r(8) = Aeee, 9 > 0, (1.1) 


where A denotes the radius associated to 0 = 0. It gives the distance of a curve point to 
origin O in terms of 0. In cartesian coordinates, the points (x(0), y(0)) of the spiral are 


given by 


Note that when a = 90° the equiangular spiral degenerates to a circle. Of course the 
animal would not be very satisfied with a circular shell, because he could not keep growing 
inside the shell. If œa is not a right angle, then a true spiral forms, which corresponds to an 
enlargement of the shell. 

This growth process keeps the shape of the shell and is called gnomonic. In geometry, the 
gnomon (a word of Greek origin, due to Aristotle, meaning “indicator”) of a given picture is 
a second picture that, added or subtracted to the former, generates a third picture similar 
to the original one. This growth pattern is so common that it is referred by many as a 
“law of nature”. The majority of animal horns and nails, corals and snails, among other 
examples, also follow, basically, equiangular spirals. 

The next figure shows three of the cases that may occur. The first example is typical of 
a cone, the second one of the nautilus shell and the third one of the clam of a bivalve. On 
the right, an horizontal section — in the case of the nautilus and the bivalve — and a vertical 
section — in the case of the cone — show the corresponding growth spiral. In each case the 
new shell material is progressively added in the aperture of the shell. 


1B quivalently, the equation may be given by log(r(9@)/A) = 0 cot a. By that reason, the equiangular spiral 


is also known as the logarithmic spiral. 
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(3) Bivalve: 


As we shall see, seashells, with their auto-similar shape, may be represented by a three- 
dimensional surface, generated by a simple equation, with some free parameters. Amazingly, 
in spite of the simplicity of that equation, it is possible to generate a great variety of seashell 
types. Which ones? All of them! (with a very few exceptions: some live and fossil species of 
Vermicularia and fossil ammonites of the class of Didymoceras.) This shows how many of the 
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forms that appear in nature are simple consequences of the application of three-dimensional 


geometry to the basic rules of growth. 


<Mathematics is very much like poetry... what makes a good 
poem — a great poem — is that there is a large amount of 
thought expressed in very few words. In this sense formulas 
like e° +1=0 or JS; e dx = yT are poems.> 

— LIPMAN BERS 


2. THE MODEL 


The surface of a shell is a three-dimensional surface that may be regarded as the result of 
the deplacement of a curve C (the generating curve, usually an ellipse) along an helico-spiral 
H (the structural curve) ([I], [4]); the width of C increases as far as it moves along H: 


n/ 
a 


Why is H an helico-spiral? Basically because the mollusk does not enlarge its shell in 
a uniform manner: it secretes shell material faster on one side than the other of the open 
edge of the shell. 

The shape of C describes the outline of the shell sections and of the shell aperture while 
H determines the global shape of the shell. C is not always an ellipse. This is the case of 
the beautiful Japanese Wonder shell, which is generated by a triangular curve: 
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The effect of C in the final shape of the shell is illustrated in Section 4 (see also [7] for 
more details). 


The first step in the construction of the model is to fix a three-dimensional cartesian 
system of coordinates XY Z and to consider the parametric equation of the helico-spiral, 


written in polar coordinates r and @: 


Shell surface 


Let a denote the equiangular angle of spiral H, 8 the enlarging angle of spiral H (that 
is, the angle between Z-axis and the line from aperture local origin to XY Z origin) and let 
A be the size of the spiral aperture (that is, the distance from main origin of aperture at 
0=0): 
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Side view: 


Top view: e Z 
X 
X 
y 
y 
y A 
a=90" a=60" a=30" B 


Seen from above the helico-spiral H looks like a logarithmic spiral. So, using (1.1), we 
assume that, for each 0 > 0, the distance r(0) of the corresponding point (x(0), y(0), z(0)) 
of H to the origin is given by r(0) = Aef ta, Then 


2(0) = —r(6) cos 8 


and 


Therefore, the points H(0) = (x(@), y(@), z(@)) of the helico-spiral satisfy equations 
z(0) = Asin 8 cos 6 e? cota 
H(0) =< (6) = Asing sin 6 ef to (2.1) 


z(0) = —Acos 8 efte, 
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The generating curve C that determines the surface of the shell is, in most cases, an 


ellipse with parameters a (semimajor axis) and b (semiminor axis), 


L 
NLY 


that is, a curve with parametric equation 
ab cos s\ 2 sin s\2]—3 
re(s) = == |( ) +( )] ,O<s<2n. (2.2) 
[(a sin s)? + (bcos s)?]2 a b 
The width of C increases as far as it moves along H. We assume that its increasing rate 
ri(0) is the same as the one of H, that is, r;(0) = efte, Then, the equation of each C in 


polar coordinates (centered at the corresponding H(6)) is given by 


R.(6, 8) = rels) ri(0) = rels) 2, 0 > 0, O< 5 < 2r. 


Equivalently, in terms of cartesian coordinates: 


Z 
À 
x° (0, s) = cos s cos O re( s) e? + ; 
C(O, s) = y (0, s) = COS ssin 0 re(s) e? cot a 0.3) 
z< (0, s) = sin s re(s) ef cota 
y 


Finally, to obtain the equation of the shell G, it suffices to put the equations (2.3) of C 
in the corresponding points of H (ruled by equations (2.1)): 


xrS(0, s) = (A sin 8 cos 0 + cos s cos 6 re(s)) efta 
6(0,s) = H(0)+C(0,s)= < yS(0,s) = (A sin sin 0 + cos ssin 0 re(s)) et (2.4) 
z©(0,s) = (—A cos 8 + sin s re(s)) e2 t9, 


We can now easily extend this model to the more general situation where the generating 
ellipse C rotates in space. For this we specify three angles ¢,2, u that establish the orien- 
tation of C in space. They measure the rotation of the curve around, respectively, a vector 
orthogonal to the plane of the ellipse, the OZ axis and its horizontal axis, as the following 
picture illustrates: 
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Z 


Generating curve 


Helico- Spiral 


Shell surface 


e Case 1 (rotation ¢): it suffices to replace s by s + ¢ in each sin and cos function in 


(2.3). 


e Case 2 (rotation Q): it suffices to replace 6 by 0 +Q in each sin and cos function in 


(2.3). 


) 


SN 


Thus, equations (2.3) are replaced by 
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x (0, s) =cos(s + ¢) cos(@ + 2) re(s) efta 
C2(8,s) = < y§(0,s) = cos(s + $) sin(@ + 2) re(s) e8 (2.5) 
26 (0,8) = sin(s + $) re(s) e? 0t 2 
that is, the equations of of the shell are now given by 
x$ (0, s) = (A sin 8 cos 6 + cos(s + ¢) cos(6 + Q) re(s)) e8 
62(6,s) = y$ (8, s) = (A sin sin 0 + cos(s + ¢) sin(@ + Q) r-(s)) ef cota (2.6) 


z9(0,8) = (—A cos B +sin(s + Q) re(s)) efta, 


e Case 3 (rotation u): Observing the ellipse C2 by profile and the result of its rotation 
of angle p, let C3(0,s) = (x$ (0, s), yf (0, s), 2§(0,8)) be the point on the new ellipse 
C3 corresponding to C2(0, s): 


C,(0, 8) 


H(9) = (x(6), y(@),2(8)) 


Then 
2 (0,8) = z(0) + 2§(0, s) cos u = (—A cos 8 + cos psin(s + o) re(s)) e°°**. (2.7) 


Further, looking from above we have 


from which we conclude that 
a$(0,8s) = 29(0,s)-—a=2$(6,s) —[z$(0,s) sin u] sin(8 + 9) (2.8) 
y5(0,s) = yy (0,8) +b = y9 (0, s) + [25 (0, s) sin u] cos(8 + 9). (2.9) 
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Finally, taking also into account the direction D of coiling (1 if dextral and -1 if sinistral), 
we get the parametric equations that describe the surface of a shell generated by an ellipse 
C, with equation re(s), moving along an helico-spiral (D. Fowler, H. Meinhardt and P. 
Prusinkiewicz [4], M. B. Cortie 


D[Asin 8 cos 6 + cos(s + ¢) cos(0 + Q) re(s) — sin p sin(s + p) sin(0 + Q) re(s)] ef 8% 


SG 
Y3 (0,5) — 
[Asin Gsin 8 + cos(s + ¢) sin(@ +2) re(s) + sin u sin(s + ¢) cos(@ + Q) re(s)] e? t2 


z5(0, 5) =[-Acos@ + cos u sin(s + 6) re(s)] €9**. 


3. THE COMPLETE MODEL: ALLOWING NODULES AND SPIKES 


In case we want to add nodules, bumps and spikes to the shell, it suffices to insert some 
additional parameters: N (the number of nodules along a complete revolution of 0), Wi 
(the length of each nodule along the generating curve), W2 (the length of each nodule along 
the helico-spiral), L (the height of each nodule) and P (the angle that indicates the position 
of the nodule in the generating curve): 


Wi 


W2 


oa 


Top view Side view 


From [2] (cf. [5]), if a surface sculpture term is superimposed in the ellipse C, it suffices 
to replace the equation re(s) of the ellipse by re(s) + rn(s,@), (0 < s < 27, 0 > 0), where 


2 Alternatively, C may be any other planar curve such as, for example, a cardioid with equation re(s) = 


a(l — cos s5). 
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rn(s,0) = 0 in case W1 = 0 or Wa = 0 or N = 0, and 
2(s—P) 21(0) 
mMehapelm tim)! y pe im ( 28), (3.1) 


otherwisd?| 


In summary, the equations of the shell depend on 14 parameters (the latter four are only 


of relevance in case L is nonzero): 


D, A, a, 8, u, Q, 6, a, b, L, P, W1, Wo, N. 


«Mathematics, rightly viewed, possesses not only truth, but 
supreme beauty — a beauty cold and austere, like that of sculp- 
ture, without appeal to any part of our weaker nature... sub- 
limely pure, and capable of a stern perfection such as only the 
greatest art can show.> 


— BERTRAND RUSSELL 


4. EXAMPLES 


The model above, incorporating (3.1), can be easily written in Mathematica in the fol- 
lowing wayi} 


1[n_] [theta_]=(2Pi/n)((n theta/2Pi)-IntegerPart[n theta/2Pi]); 
k[a_,b_,11_,p_,wi_,w2_,n_][s_,theta_] = If[w1=0| |w2=0] |n=0,0, 
11 Exp[-(2(s-p) /w1)*2-(21[n] [theta] /w2)°2]]; 
h[a_,b_,11_,p_,w1_,w2_,n_] [s_,theta_]=(1/(Sqrt[(Cos[s] /a)~2+(Sin[s]/b)~2])) 
+ k[a_,b_,11_,p_,wi_,w2_,n_][s_,theta_]; 
x[d_,aa_,beta_,phi_,omega_,mu_,alpha_,a_,b_,11_,p_,wi_,w2_,n_] [theta_,s_]= 
d(aa Sin[beta]Cos[theta]+h[a,b,11,p,wi,w2,n] [s, theta] 
(Cos [st+phi] Cos [thetatomega] -Sin [mu] Sin[stphi] Sin[thetatomega] )) 
Exp[theta Cot [alpha] ]; 
y [aa_,beta_,phi_,omega_,mu_,alpha_,a_,b_,11_,p_,w1_,w2_,n_] [theta_,s_]= 
(-aa Sin[beta] Sin[theta]-hla,b,11,p,wi1,w2,n] [s,theta] 
(Cos [s+phi] Sin [thetatomega] +Sin [mu] Sin[stphi] Cos [thetatomega] )) 
Exp[theta Cot [alpha] ]; 
z[aa_,beta_,phi_,mu_,alpha_,a_,b_,11_,p_,wi_,w2_,n_] [theta_,s_]= 
(-aa Cos[beta]+h[a,b,11,p,w1,w2,n] [s,theta]Sin[s+phi] Cos [mu] ) 
Exp[theta Cot [alpha] ]; 


3Note that this model contains the previous case of shells with no nodules: it suffices to make L = 0. 
4Since the XY Z referential we use in our model has opposite orientation to the one used by Mathematica, 


we change the sign in coordinate y. 


SEASHELLS: THE PLAINNESS AND BEAUTY OF THEIR MATHEMATICAL DESCRIPTION 13 


By choosing appropriate values for the various parameters and ranging 0 from Omin to 


O@maz and s from smin to smax 


Top view Side view 


we can now plot the surface of any type of shell. 

For example, the Boat Ear Moon (cf. [8], p. 78), from the class of Gastropods, can be 
generated by parameters D = 1, A = 25, a = 83°, B = 42°, u = 10°, Q = 30°, ¢ = 70°, 
a = 12, b = 20, L = 0, P = 0, Wi = 1, Wo = 1 and N = 1: 


With [{d=1,aa=25 ,beta=42Degree ,phi=70Degree ,omega=30Degree ,mu=10Degree, 
alpha=83Degree ,a=12,b=20,11=0 ,p=0 ,w1=0 ,w2=0 ,n=0}, 
ParametricPlot3D[{x[d,aa,beta,phi,omega,mu,alpha,a,b,11,p,w1,w2,n] [theta,s], 
y[aa,beta,phi,omega,mu,alpha,a,b,1l1,p,w1,w2,n] [theta,s], 
z[aa,beta,phi,mu,alpha,a,b,1l1,p,w1,w2,n] [theta,s]}, 
{theta, -4Pi, 4Pi}, {s, -270 Degree, 90 Degree}, 
Boxed -> False, Axes -> False, PlotPoints -> {100,20}, 
PlotRange -> All, ViewPoint -> {-1,-3,0.5}]]} 
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We conclude with a few more examples of shells that may be found in [3] and that can 
be modelled this way (in [7] you may manipulate the figures, with more detail). The most 
surprising example is the one of bivalves, with their umbones (protuberances in the top) 


wonderfully plotted! 


The Precious Wentletrap 
(D =1, A= 90, a = 86°, 8 = 10°, p = 5°, Q = 1°, ġ = —45°, a = 20, b = 20, L = 14, P = 40, 
W, = 180, W2 = 0.4, N = 180): 


Turritella, ancilla and whelk: 
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Two top shells (the Tiger Maurea and the Commercial Trochus): 


Tun shell, olive and horse conch: 


es. 


Two limpets (the Common European Limpet and the Blue-rayed Limpet): 
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Two bivalves (the European Glycymeris and the Lyonsia): 


Argonauta (Cephalopod): 
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Worm shell, elephant tusk (Scaphopod) and spirula: 


<The scientist does not study nature because it is useful; he 
studies it because he delights in it, and he delights in it because 
it is beautiful. If nature were not beautiful, it would not be 
worth knowing, and if nature were not worth knowing, life 
would not be worth living.> 


— HENRI POINCARE 
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